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I.  INTRODUCTION 

Because  of  the  development  of  oil  resources  in  the  Arctic,  it  has 
become  necessary  to  devise  methods  of  dealing  with  accidental  spills  of 
oil  beneath  and  near  sea  ice.  One  of  the  first  steps  involved  in  clean¬ 
ing  up  an  under-ice  oil  spill  is  determining  its  exact  location.  If  a 
means  of  rapidly  locating  an  under-ice  oil  spill  can  be  developed, 
appropriate  recovery  techniques  that  capitalize  on  that  information  can 
then  be  considered.  At  the  request  of  the  U.S.  Coast  Guard,  the  Applied 
Physics  Laboratory  undertook  a  study  to  determine  the  feasibility  of 
using  underwater  acoustics  to  locate  under-ice  oil  spills.  That  study, 
which  is  summarized  here,  concentrated  on  the  problem  of  localizing  an 
under-ice  oil  spill  whose  general  position  is  known.  First  we  measured 
the  surface  backscattering  strength  of  sea  ice  at  100-300  kHz  at  low 
grazing  angles.  We  then  used  the  results  to  calculate  the  distance  at 
which  oil  spills  should  be  detectable. 


I I .  SUMMARY 

At  low  grazing  angles,  acoustic  backscatter  from  the  relatively 
rough  ice-water  interface  should  be  substantially  larger  than  that  from 
the  smooth  oil-water  interface.  To  obtain  quantitative  data  on  this 
difference,  we  conducted  experimental  investigations,  in  both  the  labora¬ 
tory  and  the  field,  to  determine  the  backscattering  strength  of  flat  sea 
ice  in  the  frequency  range  of  f  -  100-300  kHz,  and  for  grazing  angles  of 
0  =  10°-45° . 

The  surface  backscattering  strength,  S$,  was  found  to  be 
Ss  =  -72  +  25  log  f  +  15  log  sin  9,  dB, 
with  a  standard  deviation  of  5.25  dB. 

The  backscatter  from  an  oil  layer  under  ice  was  too  small  to  be 
measured  at  small  grazing  angles  with  our  system. 

Using  the  empirical  formula  of  Eq.  (1),*  a  sample  calculation  was 
made  which  showed  that  an  under-ice  sonar  could  effectively  map  the 
location  of  spilled  oil  under  ice  in  an  area  several  hundred  meters  in 
diameter.  The  required  resolution,  range,  and  display  tradeoffs  have 
yet  to  be  calculated,  but  the  necessary  acoustic  design  information  has 
been  obtained. 


* 

See  p.  4. 
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III.  BACKGROUND 

Sea  ice  is  a  very  inhomogeneous  material.  Freezing  starts  at  the 
sea  surface  with  small  disk-shaped  particles  (platelets)  which  grow  in  a 
snowflake- like  pattern;  the  perpendicular  line  through  the  center  of  the 
disk  is  the  "C"  axis  of  the  ice  crystal.  These  particles  initially 
float  with  their  C  axis  vertical.  As  freezing  continues  and  the  crystals 
begin  to  crowd  each  other,  they  eventually  tip  so  that  the  C  axis  becomes 
horizontal.  In  the  process,  brine  is  concentrated  between  the  platelets. 
The  freezing  process  continues  with  the  downward  growth  of  the  ice 
platelets,  the  eventual  freezing  up  of  the  space  between  the  platelets, 
and  the  displacement  and  some  entrapment  of  the  brine  concentrate.  We 
have  developed  a  casting  technique  that  allows  us  to  make  accurate 
impressions  of  the  ice  surface,  or  "skeletal  layer."  This  layer  is 
typically  2  or  3  cm  deep  and  is  quite  fragile.  Figure  1  is  a  photograph 
of  a  casting  taken  from  a  sample  under-ice  skeletal  layer.  The  platelets 
are  about  0.8  mm  apart. 


Figure  1.  A  easting  of  the  under-iee  skeletal  layer. 
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The  dispersion  characteristics  of  oil  that  is  spilled  under  sea  ice 
have  been  extensively  studied.  Experiments  by  the  U.S.  Coast  Guard1*2 
in  the  early  '70s  provided  an  excellent  background  for  the  more  defini¬ 
tive  experiments  and  studies  that  were  sponsored  by  the  Department  of 
the  Environment,  Canada,  in  conjunction  with  the  Beaufort  Sea  Project.3'6 
An  excellent  summary  of  this  knowledge  was  given  by  Lewis7  in  1976. 

Oil  that  collects  at  the  ice-water  interface  is  always  less  dense 
than  seawater.  If  the  underside  of  the  ice  is  perfectly  flat,  the  oil 
disperses  horizontally  under  gravitational  force  until  an  equilibrium 
thickness  of  slightly  less  than  a  centimeter  is  reached.  The  ice-water 
interface,  however,  is  generally  not  flat.  Small-scale  relief  in  the 
under-ice  surface  is  usually  related  to  the  insulating  effects  of  the 
varying  depths  of  snow  cover.  In  the  spring,  the  amplitude  of  the  re¬ 
sulting  undulations  can  easily  be  20  to  30  cm  with  wavelengths  of  5  to 
10  m,  corresponding  to  snow  drift  characteristics.  Oil  will  thus  tend 
to  collect  in  pools. 

Consider  the  case  of  oil  spilled  under  ice  during  the  freezing  sea¬ 
son.  As  the  oil  lies  against  the  ice  it  impedes  the  flow  of  heat  through 
the  ice  from  the  water  to  the  air.  Because  oil  has  a  thermal  conductiv¬ 
ity  about  one-fifteenth  that  of  sea  ice,  its  presence  reduces  the  temper¬ 
ature  gradient  previously  existing  in  the  ice,  as  shown  in  Figure  2. 


Figure  2.  Effect  of  an  oil  layer  beneath  the  ioe  on  the  freezing 
proae88. 
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Because  the  rate  of  heat  transfer  is  fastest  at  the  edge  of  the  oil 
pool,  a  lip  of  ice  builds  up  which  further  restrains  the  horizontal 
movement  of  the  oil.  Eventually,  a  new  layer  of  ice  begins  to  form  on 
the  lower  surface  of  the  oil,  and  the  oil  becomes  encapsulated  in  the 
sea  ice.  The  time  required  for  encapsulation  depends  on  ice  thickness, 
air  temperature,  and  oil  pool  thickness.  In  the  NORCO  tests,3  this 
period  varied  from  about  5  days  in  the  fall  to  10  days  in  the  spring. 

The  time  limit  for  localizing  and  recovering  the  oil  before  it  is  en¬ 
capsulated  is  therefore  short  during  the  freezing  season.  (During  the 
melt  season,  the  oil  tends  to  rise  through  the  brine  channels  created  in 
the  ice  by  the  freezing  process,  resulting  in  collection  problems  different 
from  those  in  other  seasons.) 

The  idealized,  or  assumed,  picture  of  an  under-ice  oil  spill  there¬ 
fore  is  a  patch  of  oil  with  a  smooth,  flat  surface  surrounded  by  ice  of 
relatively  rough  texture.  The  feasibility  of  locating  the  oil  acous¬ 
tically  depends  on  there  being  significant  backscatter  from  the  s^a  ice 
around  the  oil  patch,  in  order  to  define  its  boundaries.  This  is  analo¬ 
gous  to  the  usual  radar  piloting  problem  of  identifying  navigable  water 
(the  area  with  no  return  signals)  as  opposed  to  the  shoreline,  buoys, 
etc.  (the  areas  with  strong  return  signals). 

The  localization  system  would  consist  of  an  acoustic  transducer 
below  the  ice-water  interface,  which  would  directionally  transmit  and 
receive  pulses  in  such  a  way  as  to  create  a  plan  view  map  of  the  local 
area.  The  map,  analogous  to  a  radar  display,  would  show  darkness  in 
regions  of  low  reflectivity,  presumably  the  oil  pool,  whereas  areas  with 
high-level  returns  (brightness)  would  indicate  regions  of  sea  ice. 


IV.  CALCULATION  OF  SURFACE  BACKSCATTERING  STRENGTH 

The  parameter  that  characterizes  the  backscattering  properties  of 
a  surface  is  called  the  surface  scattering  strength,  Ss.  Using  Urick's8 
notation,  Ss  is  related  to  other  variables  by 

RL  =  SL  -  40  log  r  -  2ar  +  Ss  +  10  log  A,  dB//uPa,  (1) 

where 

RL  is  the  received  level  (dB//yPa) 

=  20  log  V.  -  S  .  where 
r  x 

V  =  received  voltage 

*  receiving  sensitivity  (dB//V/yPa) 

SL  is  the  source  level  of  the  acoustic  projector  (dB//yPa  8  1m) 
=  20  log  Vt  ♦  Tv,  where 

Vt  ■  transmitted  voltage 

Tv  =  transmitting  sensitivity  (dB//pPa/V  91m) 
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r  is  the  range  (m) 


a  is  the  absorption  coefficient  (dB/m) 


A  is  the  equivalent  area  from  which  the  signal  is  returned 
during  a  short  time  interval  (m2) 

a  CT 

r  2  cos  0 

for  short  acoustic  pulses  and  small  grazing  angles, 
where 


(2) 


$  =  equivalent  ideal  full  horizontal  beam  width  (radians) 
for  uniform  two-way  response8 

c  =  velocity  of  sound  in  water  (m/s) 

x  =  pulse  length  (s) 

6  =  grazing  angle  (deg)  . 


Equation  (2)  can  also  be  rewritten  as 

10  log  A  =  10  log  r  +  10  log  $  +  10  log  -  ^  .  (3) 

Substituting  (3)  into  (1),  we  obtain 

RL  =  SL  -30  log  r  -  2ar  +  S«.  +  10  log  3>  +  10  log  _-CT — .  (4) 

5  2  cos  0 

We  then  solve  (4)  for  Ss : 

Ss  =  RL  -  SL  +  30  log  r  +  2ar  -  10  log  $  -  10  log  2-^g  .  (5) 

The  surface  scattering  strength  is  determined  experimentally  from 
measurements  of  the  quantities  in  Eq.  (5). 


V.  EXPERIMENTAL  METHODS 

A  transducer  was  designed  to  measure  backscatter  from  the  under-ice 
surface  over  a  wide  frequency  range.  It  consisted  of  a  22  x  22  element 
array  cut  from  an  8.26  cm  square  ceramic  plate  (Figure  3).  The  elements 
of  the  center  4x4  array  were  electrically  connected  to  a  common  lead, 
as  were  those  of  each  successive  row  grouped  around  the  center  array, 
resulting  in  10  electrical  leads  (plus  a  common  back-plane  connection) 
which  could  be  combined  in  various  ways  to  obtain  different  transducer 
operating  characteristics.  The  assembled  transducer,  potted  in  polyure¬ 
thane  and  mounted  on  a  test  fixture,  is  shown  in  Figure  4. 

The  transducer  could  be  operated  in  either  of  two  modes: 
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Figure  3.  Array  transducer  of  22  x  22  elements  used  for  backscattering 
measurements . 


Figure  4. 

Transducer  assembled  and  mounted  on 
the  movable  am.  The  electrical 
leads  are  brought  out  through  a 
rubber  tube. 
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(1)  Starting  with  the  4x4  array  in  the  center,  more  and 
more  squares  of  elements  could  be  summed  to  form  larger 
arrays.  For  each  array,  a  frequency  could  be  selected 
such  that  the  beam  had  a  particular  width  so  that  the 
beamwidth  could  be  held  constant  regardless  of  which 
array  was  used.  In  this  mode,  however,  the  beams  had 
relatively  large  side  lobes. 

(2)  All  22  x  22  elements  could  be  used,  with  proper  shading, 
to  produce  beam  patterns  with  suppressed  side  lobes.  The 
width  of  the  beam  was,  however,  different  at  different 
frequencies . 

Mode  2  was  used  because  side  lobe  suppression  was  highly  desirable,  con¬ 
sidering  the  confined  space  for  our  laboratory  experiments.  The  beam 
width  as  a  function  of  frequency  for  this  mode  is  shown  in  Figure  5. 

For  reasons  of  azimuthal  resolution,  relatively  narrow  beams  should  be 
used.  Therefore  only  frequencies  of  100,  145,  200,  and  300  kHz  were 
employed.  The  beam  patterns  at  th«se  frequencies  are  shown  in  Figure  6. 


Figure  5.  Variation  in  beam  width  of  shaded  22  x  22  element  transducer . 
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In  the  laboratory,  testing  took  place  in  an  81  cm  wide  *  180  cm 
long  x  117  cm  deep  insulated  tank  in  a  temperature-controlled  cold  room. 
The  tank  was  filled  with  salt  water  at  32 %e  salinity.  The  transducer 
was  mounted  on  a  circular  track  which  was  positioned  so  that  the  center 
of  the  acoustic  beam  intercepted  the  ice-water  interface  at  the  same 
point,  regardless  of  grazing  angle.  The  entrance  hole  was  sawed  through 
the  ice  after  the  desired  ice  thickness  was  reached.  A  schematic  of  the 
experimental  setup  is  shown  in  Figure  7,  and  a  photograph  of  the  tank  in 
the  cold  room  is  shown  in  Figure  8. 

The  ice  was  frozen  at  different  rates  to  simulate  the  freezing 
rates  that  would  be  experienced  during  different  arctic  seasons.  The 
formation  of  a  skeletal  layer  was  observed  in  all  cases.  For  each  case, 
backscatter  from  the  ice  was  measured  at  grazing  angles  of  10°,  15°, 

20°,  30°,  and  45°.  At  each  angle,  80  us  long  pulses  were  transmitted  at 
100,  145,  200,  and  300  kHz.  The  backscattered  signals  were  acquired  and 
stored  by  a  Nicolet  digital  oscilloscope  on  floppy  diskettes,  and  later 
recalled  for  analysis. 

The  measurement  process  was  repeated  on  an  arctic  field  trip  to 
the  Beaufort  Sea  in  October  and  November  1982.  The  same  transducer  was 
used  as  before  but  the  grazing  angle  was  adjusted  with  a  parallelogram 
linkage.  The  transducer  was  3  m  below  water  level  and  looked  upward  at 
60  cm  thick  ice.  The  ensonified  area  was  substantially  larger  than  that 
in  the  laboratory  tests.  Measurements  were  made  in  two  arbitrarily 
selected  directions,  about  60°  apart  in  azimuth. 


COLD  ROOM 


Figure  7 .  Block  diagram  of  experimental  setup  for  laboratory  measurement 
of  ice  backscatter. 
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j  VI .  RESULTS 

I 

j  Surface  scattering  strengths,  Ss,  were  computed  using  Eq.  (5).  The 

results  are  shown  in  Table  I.  The  mean  and  standard  deviation  of  Ss  are 
j  given  in  Table  II. 

|  It  has  been  found  that  Lambert’s  law 

Ss  =  K  +  10  log  sin20 

in  optics  is  also  a  good  model  for  describing  backscattering  of  sound 
from  rough  sea  bottoms.8  Because  the  skeletal  layer  is  somewhat  analo¬ 
gous  to  a  rough  sea  bottom,  we  used  this  model  to  formulate  the  following 
equation 

Ss  =  S0  ♦  B  log  f  +  C  log  sin  0,  f  in  kHz,  (6) 

i 

■  i 
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Table  I.  Measured  under-ice  surface  backs  oat tering  strength 
in  decibels. 
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-16.0 

-17.5 

-8.7 

-47.8 

-39.5 

-51.8 

-25 . 5 

-55.1 

-3.3.9 

-28.2 

-20.2 

-29.9 

-25.5 

-20. o 

-18.2 

-29.0 

-23.0 

-24.4 

-15.2 

-2b. 0 

-24.8 

-18.4 

-31.1 

-22.9 

-25.1 

-15.5 

-21.0 

-21.7 

-14.0 

-13.2 

-23.8 

-13.9 

-20.9 

-15.5 

-22.5 

-  lo .  8 

-  lo .  5 

-14.6 

-19.1 

-15.5 

-18.5 

-21.8 

-21.7 

-21.8 

-20.7 

-22.8 

-21.0 

-16.9 

-17.0 

-27.4 

-22.8 

-20.5 

-14.2 

-22.5 

-21.2 

-16.2 

-19.2 

-22.9 

-19.1 

-12.1 

-26.8 

-26.7 

-18.1 

-17.9 

-27.7 

-18.0 

-14.7 

-16.9 

-22.5 

-18.0 

-15.7 

-13.2 

-18.5 

-25.2 

-16.8 

-12.5 

-22.7 

-18.9 

-14.9 

-10.8 

--  indicates  no  data. 

Laboratory  set  1  was  made  with  the  transducer  in  mode  1,  and  was 
not  used. 
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Table  II.  Mean  and  standard  deviation  of  Sg  values  in  Table  I. 


Frequency  (kHz) 


Grazing 

Angle 

e 

100 

145 

200 

300 

Ss 

0 

Ss 

a 

Ss 

0 

Ss 

o 

10 

-35.27 

6.67 

-31.29 

-26.15 

-21.27 

3.94 

15 

-30.39 

6.57 

-27.36 

7.54 

-22.69 

7.68 

3.62 

20 

-28.22 

4.61 

-22.67 

5.23 

5.30 

-16.64 

30 

-24.83 

5.92 

-22.57 

4.76 

3.91 

2.52 

45 

-25.30 

4.48 

-20.21 

3.53 

-18.84 

2.68 

-12.96 

3.81 

as  a  model  for  backscattering  from  an  under-ice  surface.  The  second 
term  has  been  added  to  account  for  an  observed  frequency  dependence. 
Variables  S0,  B,  and  C  are  to  be  determined  from  experimental  data. 

Applying  least-squares  regression  to  the  average  Ss  in  Table  II  as 
a  function  of  log  f  and  log  sin  6,  we  obtained  -72,  25,  and  IS  for  S0, 

B,  and  C,  respectively.  The  standard  deviation  of  all  the  data  in 
Table  I  from  the  Ss  given  by  Eq.  (6)  was  5.25  dB. 

The  means  of  all  data  taken  at  106,  15°,  20°,  30°,  and  45°  are 
shown  in  Figure  9  along  with  plots  of  Eq.  (6).  It  is  clear  that  the 
equation  is  a  good  representation  of  the  data. 

A  sample  of  Prudhoe  Bay  crude  oil  was  provided  by  the  ARCO  refinery 
at  Ferndale,  Washington.  The  sample  had  a  density  of  0.87  g  cm"'*.  Its 
sound  velocity  of  1500  m  s-1,  at  seawater  freezing  temperatures,  was  4% 
to  5%  greater  than  that  of  the  seawater.  According  to  Snell's  law, 
this  indicates  a  critical  grazing  angle  of  ~17°.  At  any  smaller  angle, 
the  acoustic  energy  incident  on  an  oil  surface  would  theoretically  be 
reflected  fully  forward,  and  no  energy  would  enter  the  oil  layer  to 
produce  scattering  from  the  ice  above.  In  our  laboratory  tests,  these 
effects  were  observed.  At  grazing  angles  greater  than  critical,  the 
observed  backscatter  from  the  ice,  through  the  oil  layer,  was  about  the 
same  as  when  the  oil  was  absent.  The  backscatter  at  small  grazing 
angles  depends  mainly  on  the  roughness  of  the  surface  and  is  much  less 
from  the  smooth  oil-water  interface  than  from  the  rough  ice-water  interface. 
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Fig'xre  9.  Average  of  measured  Se  vs  empirical  equation 
S8  *  -72  +  25  log  f  +  15  log  sin  6. 
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VII.  SAMPLE  CALCULATIONS 

Because  acoustic  differentiation  between  ice  and  oil  appeared 
possible  at  low  grazing  angles,  we  next  explored  the  question  of  whether 
the  scattering  strength  of  the  ice  surface  is  large  enough,  compared 
with  ambient  or  system  noise,  to  yield  adequate  survey  ranges. 

For  this  investigation,  we  used  a  7.2  cm  diameter  circular  piston 
transducer.  The  transducer  characteristics  and  ambient  conditions  were 
as  follows: 

f,  operating  frequency,  200  kHz 

a,  absorption,9*10  0.0387  dB/m 

beam  width  between  -3  dB  points,  6° 

directivity  index,  30  dB 

efficiency,  50%,  or  -3  dB 

4>,  equivalent  two-way  beamwidth,  0.157  rad 

SL,  source  level,11  212.6  dB//yPa  at  1  m 

Sx,  receiving  sensitivity,  -190  dBV/yPa 

BW,  receiving  bandwidth,  20  kHz 

t,  pulse  length,  100  ys 

DT,  detection  threshold,  10  dB 

Ne,  electronic  noise,  6  nV/>4fz 

Na,  ambient  noise,12  -2  dB//yPa  in  1  Hz  band. 

For  ambient  and  electronic  noise,  the  total  equivalent  in-band  acoustic 
noise  level,  NL,  is 

Ambient  NL  *  Na  +  10  log  BW  -  41  dB  , 

Electronic  NL  *  20  log  Ne  ♦  10  log  BW  -  Sx  *  68.6  dB  . 

Therefore,  the  limiting  noise  contribution  is  from  the  electronics  and 
we  used  the  latter  value  for  NL. 
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The  applicable  sonar  equation  in  this  case  is 
SL  -  30  log  r  -  2ar  +  Ss  +  10  log  $  +  10  log  y  ^  Q  -  NL  =  DT.  (7) 

In  this  equation,  r  represents  the  slant  range  and  therefore 


where  D  is  the  depth  of  the  sonar  transducer. 

Substituting  this  relationship  into  Eq.  (6),  we  obtained: 

Ss  =  SQ  ♦  B  log  f  ♦  C  log  |  .  (8) 

Substituting  Eq.  (8)  for  Ss  in  (7),  simplifying  and  rearranging 
terms,  and  noting  cos  0=1  for  small  0,  we  obtained 

SL  +  B  log  f  +  C  log  D  +  10  log  <2>  +  10  log  (ct/2)  -  NL  -  DT  +  S0 

=  (30  +  C)  log  r  +  2ar  .  (9) 

We  consider  this  equation,  with  SQ  =  -72  dB,  to  represent  the  50% 
confidence  level;  for  the  84%  confidence  level,  the  first  term  in  the 
equation  would  be  reduced  by  1  o  (to  -77.25  dB)  and  for  the  97.75% 
confidence  level,  by  2  o  (to  -82.5  dB) ,  etc. 

Applying  regression  values  of  B  and  C  and  the  sample  parameters 
given  previously  to  Eq.  (9)  and  solving  for  r,  we  obtained  the  following 
effective  ranges  for  a  transducer  depth  of  3  m  at  the  specified  confidence 
levels : 


Confidence 

Level 

C%2 


50 

84 

98 


Detection 

Range 

(m) 


149.3 

125.1 

103.6 


These  are  quite  usable  distances.  Within  these  ranges,  for  grazing 
angles  less  than  critical,  backscatter  from  the  ice  should  show  up  well 
compared  with  that  from  the  oil.  Note  that  for  these  distances  we 
are  in  fact  extrapolating  Ss  to  angles  less  than  10°.  To  compensate  for 
the  uncertainty  resulting  from  extrapolation,  a  high  nominal  confidence 
level  can  be  required  if  necessary. 
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In  general,  the  design  of  an  operational  survey  system  would  be 
relatively  straightforward,  once  the  parameters  of  the  system  are  defined. 
The  type  of  display  system  selected  will  influence  the  design  of  the 
beam  pattern.  A  very  narrow  horizontal  beam  is  desirable  for  good 
resolution  at  long  range.  However,  this  reduces  the  ensonified  area  and 
thus  the  backscattered  energy.  Similarly,  a  higher  operating  frequency 
will  produce  a  higher  Ss  as  well  as  higher  sound  absorption,  each  affect¬ 
ing  the  detection  distance  in  opposite  ways.  Optimization  studies  of 
parameters  such  as  these  would  be  the  next  step  toward  the  construction 
of  an  acoustic  system  for  locating  under-ice  oil. 
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